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ABSTRACT 


Asymptotic  methods  are  used  to  find  approximate  sol¬ 
utions  of  the  acoustic  wave  equation  in  a  medium  where 
the  velocity  is  a  continuously  variable  function  of  one 
coordinate.  It  is  shown  that,  when  the  velocity  func¬ 
tion  has  a  minimum,  undamped  normal  mode  solutions 
exist,  and  that  such  solutions  are  closely  analogous  to 
the  internally  reflected  waves  in  the  case  of  a  medium 
made  up  of  discrete  layers.  By  converting  the  sum  of 
the  high-  order  normal  modes  into  an  equivalent  inte¬ 
gral,  it  is  shown  that  superposition  of  these  modes 
leads  to  geometrical  ray  theory  modified  by  diffraction 
in  a  manner  that  Jiay  be  computed  from  the  incomplete 
Fresnel  and  Airy  integrals. 
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DIFFRACTION  EFFECTS  IN  THE  PROPAGATION 
OF  COM PRES SI ON AL  WAVES  IN  THE  ATMOSPHERE* 


I  NTRODIICTION 


Following  F.  J.  W.  Whipple’s  explanation  of  the  zones  of  silence  and  of  abnormal 
audibility  of  sounds  from  large  explosions  as  an  effect  produced  by  the  temperature  struc¬ 
ture  of  the  atmosphere,  a  number  of  estimates  of  the  temperature  versus  altitude  function 
at  altitudes  in  the  range  from  25  to  60  km  hare  been  made2  on  the  basis  of  travel  time 
measurements.  As  in  the  analogous  problem  in  seismology,  the  determination  of  the  veloc¬ 
ity-altitude  function  from  the  observed  time-distance  curve  has  been  based  on  ray  geometry. 
The  ray  concept  is,  of  course,  valid  and  useful  if  the  frequencies  involved  are  sufficient¬ 
ly  hig^i,  but  it  does  not,  in  itself,  provide  any  criterion  for  determining  how  high  the 
frequency  must  be.  In  this  report  we  shall  consider  diffraction  effects  by  applying  a 
method  suggested  by  Pekeris3  to  find  an  approximate  solution  of  the  wave  equation  for  the 
case  in  which  the  velocity  is  a  function  of  altitude  only  and  has  a  single  minimum.  Al¬ 
though  the  discussion  will  be  carried  through  with  particular  reference  to  propagation  in 
the  atmosphere,  it  will  be  apparent  that  the  method  is  equally  applicable  to  the  long  range 
sound  transmission  in  the  ocean  which  is  associated  with  the  existence  of  a  velocity  mini- 
uusn  ("Sofar"  channel)  and,  with  some  modification,  to  the  propagation  of  seismic  waves  in 
cases  where  a  low  velocity  layer  is  present. 


•Manuscript  racalvad  10  January  19S0 


lf.  J.  W.  Whipvla.  N  •  t  u  r  • ,  »ol  HI.  P  187;  19  23. 

2fl.  Ou  t  anbt  r  | ,  Htndbuch  dtr  G^ophyik,  »ol  9.  p  89l  193  2-  A  (unwary  of  the  earlier  work  done  on 

this  subject  is  |ieen  here. 

E.  E.  Cci,  J.  Acou  *  t .  Soc.  Am.;  rol  19,  P  832,  1947;  »ol  21.  P  6.  1949.  Recent  results  Ires  the 

Arco.  Idaho  and  Helfoland  blasts  are  ditcuased  here. 

3C.  L.  Pekeris,  "Theory  of  Propafation  of  Sound  in  a  Half  Space  of  Variable  Sound  Velocity," 
/.  icons  I  .  Soc.  Am.,  vol  18.  P  295:  1946. 


8 


1.  WAVE  EQUATION 


When  the  effect  of  gravity  and  the  consequent  variation  of  initial  density  are  not 
ignored,  the  hydrodynamic  equations  of  motion  generally  do  not  possess  irrotational  solu¬ 
tions  that  can  be  expressed  in  terms  of  a  velocity  potential.4  If  the  relative  change  of 
the  initial  density  within  a  wavelength  is  small,  however,  the  motion  may  be  approximated 
locally  by  a  velocity  potential  4>  satisfying  the  simple  wave  equation, 

c2(i  )Vad>  -  Bay/Bta  =  0  ,  (  1-01) 

with  the  particle  velocity  vector  V  and  the  pressure  perturbation  p  given  by 


V  -  -  grad  <t>  ( 1  '02 ) 

D  =  p(z)  — —  ,  (1*03) 

B  t 

where  c(z)  is  the  local  velocity  of  sotxid  and  p(z)  is  the  local  density,  both  functions 
of  the  altitude  z.  If  these  equations  are  applied  to  a  case  in  which  the  range  of  alti¬ 
tudes  involved  is  large,  it  is  found  that  the  conservation  of  energy  is  violated.  The 
acoustic  energy  flux  density  if. 

Btf 

J  =  pV  =  -p  (i)  --  grad  4>  .  (1*04) 

d  t 


Since  is  a  solution  of  Eq.  (1*01),  the  factor  -  (B<h/Bt)grad  <£  in  Eq.  (1*04)  ex¬ 
presses  the  effect  of  geos>etrical  ray  divergence  or  convergence;  and  the  factor  p(z) 
then  leads  to  a  loss  of  energy  by  a  wave  travelling  in  the  direction  of  decreasing  density, 
and  a  gain  of  energy  by  a  wave  travelling  in  the  direction  of  increasing  density.  Thus, 
even  though  the  terms  depending  on  gravity  and  initial  density  variation  are  small,  their 
cumulative  effect  may  be  very  large  when  the  range  of  altitudes  considered  is  great.  We 
may  restore  the  conservation  of  energy,  however,  while  retaining  the  simplicity  of  Eq. 
(1*01),  by  multiplying  the  right-hand  sides  of  Eqs.  (1*02)  and  (1*03)  by  a  dimensionless 
quantity  f(z)  wh*ch  is  a  slowly  varying  function  of  z  in  the  sense  that  the  relative 
change  in  f(z)  within  a  wavelength  is  small.  The  equations  of  motion  and  continuity  then 
are  satisfied  locally  to  the  same  order  of  approximation  as  that  involved  in  deriving  Eqs. 
(1*01),  (1*02),  and  (1*03).  If  we  set  f(z)  =  [p(h)/p( z)] l/ 3  where  p(h)  is  the  densi¬ 

ty  at  the  altitude,  h,  of  the  source.  Eqs.  (1*02),  (1*03),  and  (1*04)  become 


4H.  Ltwb.  Hydredynaml  cm  6th  td.;  CaabrJidc*  Unlrtrtity  Pratt:  pp  54  1-  55*. 


V 


P 


POO* 

fi(z) 


1/2 

grad  4 


[p(h)p(z)]  1/2 


B<£ 

b7 


(  1  -OS) 


( 1 -Oh) 


_  3  4 

J  -  -  p(n)  — ■ —  grad  4 

B  t 


(1*07) 


Since  J,  as  given  by  Eq.  (1*07),  no  longer  depends  on  the  altitude  except  through  the  ef¬ 
fect  of  convergence  or  divergence  expressed  by  the  function  4,  the  energy  is  conserved  at. 
all  altitudes;  also,  it  is  evident  from  the  definition  that  if  p(z)  is  a  slowly  varying 
function  of  altitude,  f(z)  is  also  a  slowly  varying  function.  The  condition  for  the  va 
lidity  of  Eqs.  (1*01),  (1*05),  and  (1*06)  then  may  be  expressed  generally  by  the  inequality 


X  d( log  p) 
2  rr  dr 


«  1 


(  1  *08 


where  X  is  the  /avelength.  The  average  value  of  d  log  p/dr  to  an  altjtudt  of  r/D  km  in 
the  atmosphere  is  about  “0.14  km*1  so  that  Eq.  (1*08)  is  equivalent  to  X<<45  km,  which 
implies  that  the  frequency  must  be  v>  »  0.0067  cycle  per  second.  The  approximate  methods 

that  we  shall  apply  to  the  solution  of  Eq.  (1*01),  however,  will  provide  a  somewhat  more 
severe  low-frequency  limitation  than  this. 


2.  NORMAL  MODE  SOLUTIONS  OF  THE  WAVE  EQUATION 


Periodic  solutions  of  Eq.  (1*01)  in  a  cylindrical  coordinate  system  with  axis  passing 
through  the  source  may  be  written  in  the  form 

4(.  r,  z,  t,k)  =  e  iwt  JQ(kr  )F(  r ,  k)  ,  (2*01) 

where  F(z,k)  satisfies  the  equation 


d  *F  cu2 

-  +  F[— - - 

dr2  c2(z) 


k2] 


0 


(2*02) 


The  boundary  conditions  are 

(1)  At  the  ground  surface  (z  =  0)  the  vertical  component  of  the  particle  velocity 

vanishes . 

(2)  V(  z )  —  0  as  z  *  ® 


*  ■*  •  ** 
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(3)  In  the  neighborhood  ol  the  source  (r  -  0,  z  -  h)  the  integrated  solution 


r.i.t) 


t,k) 


dk 


must  reduce  to  the  value  appropriate  to  a  point  source,  namely, 


(2*03) 


A  exp{*<*j[t  ~  (l/c)/ra  +  (z  ~  h)3]^ 

*  ~  /r*  +  (>  -K)’  ' 

Designating  the  regions  below  and  above  the  source  by  subscripts  1  and  2,  respectively, 
condi tion  (1)  is 


(2*05) 


condition  (2)  requires  that  Fa  be  a  particular  solution  of  Eq.  (2*02)  such  that,  with 
the  factor  eicjt,  it  represents  an  upward  travelling  wave  at  large  values  of  z;  and 
condition  (3)  is  satisfied3  if  we  impose  the  conditions  at  z  -  h 

F,(h,k)  =  F2(h,k)  (2*06) 


dz 


h 


2Ak 


(2*07) 


Let  M(z,k)  and  N(z,k)  be  two  linearly  independent  solutions  of  Eq.  (2*02)  that  be¬ 
have  asymptotically  for  large  values  of  z  (with  the  factor  ela,t)  like  downward  and  up¬ 
ward  travelling  waves,  respectively.  The  boundary  conditions  are  then  satisfied  by 


Fj(z,k) 


2AkN(h,k)(N(z,k)M'(0,k) 


N'(0.k) 


-  M(z,k) 


(2*08) 


an  < 


Fa(z,k) 


2AkN(z,k)  fN.h,k)M'(0,k) 


N'(0,  k) 


-  M(h,k) 


(2*09) 


where  the  primes  denote  differentiation  with  respect  to  z,  and  b  -  MN'  -  M'N  is  a  con¬ 
stant  by  virtue  of  the  Abelian  identity  which  holds  between  any  two  linearly  independent 
solutions  of  Eq.  (2*02). 

We  now  assume  that  the  velocity  versus  altitude  function  has  the  general  character 
shown  in  Fig.  1.  This  function  and  its  first  derivative  art  taken  to  be  continuous.  There 
is  a  variety  of  evidence  for  the  existence  in  the  atmosphere  of  a  second  temperature 
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Fig.  1.  General  character  of  anuned 
velocity-altitude  function. 


minimum  at  an  altitude  of  approximately  80  km,  followed  by  a  rise  to  indefinite  heists. 
At  very  great  altitudes,  however,  attenuation,  which  is  not  allowed  for  in  Eq.  (1*09), 
becomes  increasingly  impor tan t- -and  the  more  so  the  higher  the  frequency.  Consequently, 
an  appreciable  return  of  wave  energy  to  lower  levels  by  refraction  from  altitudes  above 
the  first  maximum  is  not  to  be  expected  except  for  very  low  frequencies.  The  assumption 
of  a  monotonically  decreasing  velocity  above  the  altitude  z2,  as  indicated  in  Fig.  1, 
thus  may  be  regarded  merely  as  a  device  for  ensuring  that  any  wave  energy  that  has  trav¬ 
elled  above  the  altitude  z2  will  not  be  returned  to  lower  levels.  The  method  which  will 
be  elaborated  below  can,  in  principle,  be  applied  to  a  velocity  function  with  any  number 
of  minima  and  maxima,  but  for  the  present  we  shall  confine  our  attention  to  functions  with 
only  one  minimum  and  one  maximum. 

Pekeris3,5  has  shown  that  the  integral  (2*03)  for  can  be  reduced  by  replacing 
J0(kr)  by  [H0(1)(kr)  +  H<2)(kr)]/2  followed  by  a  transformation  of  the  path  of  integra¬ 
tion  in  the  complex  k-plane,  to  a  sum  of  residues  taken  at  the  zeros  of  N'(0,k)  plus 
certain  integrals  around  cuts  in  the  k-plane.  The  contribution  to  <t>  of  the  residues, 
i.e.,  the  "normal  mode"  solutions,  is 


4> 

•  a 


-  277iAe 


io>t 


y  kJdQ(a)(k>r)N(g,kB)N(htk>) 
„  N(0,k)[BN'(0,k  )/3k]k.k 


(  2' 10) 


where  the  kn’s  are  the  root-*  of  N'(0,kn)  =  0  end  the  constant  b  has  been  eliminated 

by  use  of  the  relation  b  =  ~N(0, kn)M' (0,  kn) .  Since  Eq.  (2*10)  is  unaltered  by  inter¬ 

changing  z  and  h,  the  expression  is  applicable  both  above  and  below  the  level  of  the 
source . 


3C-  L-  Pakarla,  "Theory  of  propagation  of  exploaiva  lound  In  ahallow  water,  Gao f  .  Soe . Am .  nemo  !  r, 
wol  27;  1948  (.Propagation  of  Sound  in  th*  Ocean). 


- • 


~«v*rr  ■*  ~r  .*  ' 
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At  large  distances  fr oai  the  source  the  radial  factor  in  Eq.  (2*10)  behaves  like 


H0’(k.r) 


1  >/l 


rrk  r 


77 


exp[-  i  (k  r  -  -j)] 


(2*11) 


In  conjunction  with  the  factor  e1"1  this  must  represent  a  progressive  wave  travelling 
in  the  direction  of  increasing  r,  so  that  the  real  part  of  ka  mu&t  be  positive  and  the 
iaiaginary  part  negative  or  sere.  Ihose  roots  that  have  a  finite  imaginary  part  (if  any 
exist)  lead  to  an  exponential  decay  factor,  so  that  for  large  values  of  r  the  signifi¬ 
cant  terms  in  Eq.  (2*10)  are  those  for  which  ka  is  real. 


3.  APPROXIMATE  EXPRESSIONS  FOR  THE  FUNCTION  N(*,k) 


We  introduce  the  paraaieter  /3(  k )  -  'o/k  and  write 


QU.k)  r 


/3  ( k )  >  c(x) 

£(k)  <  c(  z ) 


u(z,k) 


(3*01) 


(3*02) 


where  the  lower  limit  of  integration,  a,  remains  to  be  defined.  Approximate  solutions 
of  Eq.  (2*02),  asymptotically  valid  for  sufficiently  high  frequencies,  may  be  written  in 
the  form6 


F 


1/2 


C(u) 


(3*03) 


wher^  C  is  a  linear  combination  of  Bessel  functions  of  order  1/3.  The  differential  equa¬ 


tion  for  which  Eq.  (3*03)  is  an  exact  solution  is 


F*  +  (Q2  -  S’/S)  F  =  0 


(3*04) 


*B.  1.  Loaf or.  "On  tho  Cooaoctlon  Poraatao  •»<*  t  ho  Solution*  of  th»  Wot#  Equation."  Pkjm.  Wot.; 

toI  51.  9  SS* :  1*37. 
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where 


S  =  u,/ft0‘,/J 


(  3*05) 


and 


S'/S  =  1(070) J  -  ±(CT/0)  -  n<Q/u)J 


(  3 • 06  ) 


Thus  Eq.  (3*03)  is  an  approximate  solutior  of  Eq.  (2 *02)  for  values  of  z  such  that  S'/S 
is  small  compared  to  Q3.  For  a  given  value  of  (M,  k ) ,  0  and  u  are  proportional  to  to, 
so  that  S'/S  is  independent  of  u j.  Thus  it  is  always  possible  to  find  a  value  of  to 
great  enough  to  make  Q3  >>  S"/.S  except  at  points  where  Q2  -  0  or  S'/S  becomes  infi¬ 
nite.  However,  if  we  take  the  integration  constant,  a,  at  a  first-order  zero  of  Q3 ,  it 
can  be  shown  from  Eq.  (3*06)  that  S'/S  remains  finite  at  z  =  a,  so  that  this  zero  of 
0  is  not  a  singular  point  of  Eq.  (3*04).  Thus  an  expression  of  the  form  (3*03),  which  is 
a  valid  approximation  to  a  particular  solution  of  Eq.  (2*02)  for  z  >  a,  represents  the 
same  particular  solution  for  z  <  a.  If  0  has  more  than  one  zero,  say,  z  ~  a j  and 
z  =  ajf  the  range  of  validity  of  the  solution  obtained  by  setting  a  =  in  Eq.  (3*02) 

does  not  extend  beyond  z  -  a^,  since  u(a2)  is  not  zero  and  S'/S  becomes  infinite  it 
8j.  However,  another  solution  valid  on  both  sides  of  a2  is  obtained  by  setting  a  - 
in  Eq.  (3*02).  fce  then  have  two  approximate  solutions,  both  valid  in  the  range 
af  <  z  <  »2:  namely, 


1/2 


C.(u.) 


an  < 


F.  = 


1/2 


C2(u2) 


where  u]  and  u2  are  the  values  of  u  obtained  by  setting  «  3  a,  *nd  *2»  respective 
ly,  in  Eq.  (3*02),  and  Cj  and  C2  may  be  different  combinations  of  Bessel  functions  of 
order  1/3.  So  that  both  F  and  F2  may  be  asymptotic  representations  of  the  same  solu¬ 
tion  F  of  Eq.  (2*02),  it  is  evidently  necessary  that  there  be  a  connection  between  the 
coefficients  of  the  Bessel  functions  in  C,  and  Cj.  This  connection  is  readily  estab¬ 
lished  by  noting  that  the  leading  term  in  the  asymptotic  expansicn  of  C(u)  has  the  form 


C(u)  —  [2_V/2(Aelu  +  Be*‘u) 

[tzu  J 

where  A  and  B  are  constants  (in  general  complex).  F  «nd  F2  then  will  be  asymptoti¬ 
cally  equivalent  in  the  range  at  <  z  <  a2  if  their  respective  coefficients  are  connected 
by  the  relations 
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and 


where 


AJ  =  A»  eJCp  ^ui  ”  u j ^ 

Bj  =  Bj  erp[-i(Uj  -  ua>] 


(3*07) 


is  a  constant. 

For  certain  values  of  /3(k)  the  zeros  of  Q2  may  be  of  hifd»er  order  than  the  first; 
for  ex  ample,  when  /3(k)  -  c(  or  c^  (Fig.  1),  the  zeros  of  0J  at  z  -  z^  or  z^  are 
of  the  second  order.  In  this  case  the  asymptotic  solutions  must  be  expressed  in  terms  of 
Bessel  functions  of  order  1/4  in  place  of  those  of  order  1/3;  however,  these  special 
.uses  need  not  be  considered  here. 

The  quantity  S"/S  also  becomes  infinite  if  either  c'(x)  or  c*(  z)  becomes  infi¬ 
nite,  that  is,  if  there  is  a  discontinuity  in  either  the  velocity  or  the  velocity  gradient. 
In  this  esse  Eq.  (3*03)  will  not  represent  the  same  solution  of  Eq.  (2*02)  on  both  sides 
of  the  discontinuity;  that  ia,  we  must  allow  the  coefficients  of  the  Bessel  functions  to 
change  discontinuously  at  a  discontinuity  in  c  or  c',  which  corresponds  physically  to 
the  existence  of  a  reflected  wave.  We  do  not  wish  to  complicate  the  present  discussion  by 
a  consideration  of  such  effects,  and  it  is  for  this  reason  that  the  continuity  of  c  and 
c'  was  assumed  eLove. 

In  the  particular  case  under  consideration,  we  shall  use  subscripts  1,  2,  and  3  to 
distinguish  the  parts  of  the  velocity-altitude  function  in  the  ranges  0  to  z  ,  Zj  to  z2, 
and  x,  to  infinity,  respectively.  When  k )  <  c  ,  0  has  one  zero,  z  =  a3;  when 

c  <  yS(k)  <  c  ,  it  has  three  zeros,  a  ,  a  ,  and  a  ;  when  c  <  /3(k)  <  c  ,  it  has  two 
zeros,  a  and  a  ;  and  when  /3< k )  >  c  ,  it  has  no  zeros.  To  write  explicit  expressions 

4  w  4 

for  the  function  N(z,k),  we  divide  the  (/3,  z)  -plane  into  domains  (Fig.  2).  Recalling 


Fig.  2.  Division  of  (/9,z)  plane  into 
dowiaina  defined  by  the  relative 
values  of  and  c(z)  . 
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that  for  very  large  altitudes  N  is  to  represent  upward  travelling  waves,  the  appropriate 
expression  for  this  function  in  domains  A,  R,  and  C.  is 


N 


(  3) 


(  3  -  08 ) 


u 


3 


(  3-0Q) 


Since  H  j  3  is  a  miltiple  valued  function  of  u,  it  is  necessary  to  consider  the  way  in 

which  the  phase  angles  to  be  assigned  to  the  u’s  vary  between  the  different  domains.  As 

indicated  by  F.q.  (3*01),  arg  0  ~  i)  in  domains  A,  D,  and  G,  while  in  domains  B,  C,  and  E, 
arg  0  ~  ~n/ 2 .  In  the  nei  gh  1» rhood  of  a^,  Q  is  proportional  to  (z  —  a3)1/2,  so  that  u3 
is  proportional  to  ( z  -  a^)373,  or  Q3 ;  thus  arg  u3  =  3  arg  0-  Near  aa>  Q  is  proportional 
to  ( a a  -  z)  1/2  and  u2  is  proportional  to  -Q3,  so  that  arg  u3  =  3  arg  Q  +  77.  Near  a(  we 

again  have  arg  u(  =  3  arg  Q.  The  values  to  be  assigned  to  arg  u  in  the  various  domains 

then  are  given  by  the  following  scheme: 


Domain 

arg  Q 

arg  u, 

arg  u3 

arg  a3 

A 

0 

- 

- 

0 

B 

-n/2 

- 

- 

"377/2 

C 

~n/2 

- 

—n/2 

"377/2 

D 

0 

0 

77 

- 

F. 

—n/2 

-3tt/2 

- 

- 

In  domain  G,  Q  has  no  zeros  and  we  may  choose  a  arbitrarily.  In  this  case  N  is  not 
represented  by  the  same  combination  of  Bessel  functions  on  opposite  sides  of  a,  but  the 
leading  term  in  its  asymptotic  expansion  is  the  same  or.  both  sides.  For  continuity  witJh 
domain  A  we  set  a  =  z^  in  domain  G. 

(  l )  <  2  ) 

The  leading  terms  in  the  asymptotic  expansions  of  H,/3(u)  and  H1/3(u)  for  the 
values  of  arg  u  given  in  the  above  table  are  as  follows7: 
arg  u  ~  ~  77/2,  0,  or  n 


H^**(u)  ~ *  (2/7tu)1/j  exp  [i  (u  -  5tt/12)] 


(  3-10) 


7 


G  • 


N  . 


■  a  t  a  o  n  ,  B  e  a  a  c  l  Funct  Iona .  2nd  ad.; 


Caabrid..  Il-l-.raltr  Praia;  PP  7  5  -  20  :. 
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arg  u  =  3  77  /2 

Ou>  —  (2/ttu),/j  {exp[i(u  -  577/12)]  +  exp[-i(u  +  1177/12)]} 

arg  u  =  3  77  /2,  ~  77/2,  or  0 

(  a> 

Hj/3(u)  —  (2/7tu)  1/3  exp  [~i(u  -  577/12)] 

arg  u  r  77 

H | / 11 )  ~*  (2/ttu)1/j  {exp[i(u  +  1177/12)  +  exp[~i(u  ~  577/ 12)]} 

We  also  define  the  following  positive  real  quantities: 


(3*11) 


(3*12) 


(3*13) 
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By  using  the  connection  relations  (3*07)  with  the  asynptotic  expressions  (3*10)  to  (3*13), 
we  find  the  following  expressions  for  the  function  N(i,k)  in  the  various  domains: 


In  A. 


In  B, 


1/3  2K 

'  2  3 

'./ 3 

r .( 

N  = 

H  (u  )  — * 

— 

exp  < 

“1 

LoJ 

1/33 

L"0J 

1 

t 

fu  ,3 

1/2  (2), 

'  2  1 

1  1/2 

r 

N  = 

3  1 

H  (u  )  — 

exp 

ky, 

.0 

1/3  3 

jO. 

3 

5 rrf) 

kx  (i)  ~  — 

3  12 


5m 

~Y2 


} 


(3*16) 


(3*17) 


N 


In  C, 

(VT2 

loJ 

In  n, 


exP [ky j( 


'  2' 

1/2 

TO. 

exp 

"Wj(»)  +  ky  (a  )  +  — 

J  2  3  12 


(3*18) 


N  = 


fu . 


1/ 3 


(2)  r  23  1/2  2m3 

exp[ky2(a3)]H1/3(u2)  -*  2|— J  exp|ky2(a3)  +  — j 


3  ,OOS  (kx2(x)  -  - 


(3*19) 


N 


7  «p[kya(»3)l.|expj**ikx2(a1)  +  ~J~JHi/3(u  +  ®xpCikxa(*l  )^H<1/|(«1)| 


In 

E, 

VI 

o 

(21 

rO\ 

1/2 


exp 


2m  r 

ky^a3)  +  —I  |exP 


“ky,  (  *)  ~  tkXjlaj)  +  -~| 

2  cos  kX2 ' a l )  exp| ky ,( x) 


(3*20) 


771 

4 


In  G. 


N 


I  »/2 

5tt  ' 

— 

exp 

-i(u  -  — ) 

TO; 

l  12 

(3*21) 


j; 


0  dz 


With  these  expressions  for  N,  there  are  discontinuities  in  F(z,k)  with  respect  to 

variations  of  k  at  k  =  oj/Cj  and  oj/c2.  The  trans foimatioc  of  the  path  of  integration 
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in  the  integral  expression  for  4>  leads,  in  addition  to  the  contributions  from  the  resi¬ 
dues  given  by  Eq  (2*10),  to  integrals  to  and  from  oj/c  -  i00  ,  w/Cj  -  l00  ,  and 
oVc(x)  -  iao  around  the  cuts  in  the  k-plane  shown  in  Fig.  3.  There  is  no  discontinuity 


with  respect  to  k  at  the  cut  starting  at  o>/c(z)  on  the  real  axis,  so  that  the  integral 

around  this  cut  vanishes.  The  integral  around  the  cut  starting  at  co/c  j  does  not  vanish 

because  of  the  discontinuity  in  F  between  domains  B  and  E,  but  it  may  be  shown  that 
this  integral  is  of  order  r'  1  as  compared  to  the  residue  terms  which  are  of  order 

r*1/J.  This  integral  is  therefore  negligible  when  r  is  large.  Similarly,  the  integral 

around  the  cut  starting  at  co/ c ^  does  not  vanish.  It  will  appear  presently  that  this 
integral  is  to  be  interpreted  as  representing  the  effect  of  rays  that  leave  the  source  at 
angles  steep  enough  so  that  they  are  not  refracted  back  below  the  level  z^.  This 
integral,  therefore,  is  negligible  for  large  values  of  r  when  z  <  zJ(  which  is  the 
only  case  which  concerns  us,  but  it  would  not  be  negligible  if  we  wished  to  compute  values 
of  for  z  >  ij. 


4.  THE  ROOTS  OF  N* (0,k) 


We  have 


(  4-01  ) 
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where  H  stands  for  either  of  the  Bessel  functions  H(1*  or  H(a*.  For  a  given  value 
of  /3(k)  the  expression  in  brackets  in  Eq.  (4’01)  is  independent  of  co,  while  (uQ)1/a 
is  proportional  to  co;  hence  we  may  write 


'AT 

1/ 2  1 

I _ 

H 

lo 

1/3 

. 

(uQ),/a  H.a/J(u) 


(  4*02) 


to  the  same  approximation  as  that  involved  in  writing  N  in  terms  of  H  j.  Since  the 
leading  term  in  the  asymptotic  expansion  of  H_J/3  may  be  obtained  from  that  for  HJ/3 
by  multiplying  the  eia  term  by  c77 1  / 1  and  the  e‘iu  term  by  e*wl/a  for  any  value 
of  arg  u,  the  asymptotic  expressions  for  N'(0,k)  in  the  various  domains  are  derived 
readily  from  those  for  N  by  inspection. 


In  B, 


.<  2) 


N'(o,k)  =  [u3(o)  o(o)],/a  H;a;3[u3(o)]  — 

v. 

In  E, 

N'(0,k)  =  [u^O)  Q(0)]l/J  exp[kya(a3)]  |exp  |S  kx  a(  a  J ) 


r2CK0K| 


77 


i/: 


exp 


Tri'l 

k,*,0>  '  n) 


(4*03) 


4m -) 


+  «p[>kxJ(ai)]H<^j[u|(n)]j 


r2O(0h  i/a 

f  27711 

r 

r  37ti, 

RH 

ky,(a3)  +  ~r 

<  exp 

“ky  t ( 0 )  -  ikXjfaj)  +  — 

L  77  J 

o 

l 

2  cos  kXjCaj)  exp 


3rri' 

ky ,  ( o )  -  — 

4 


} 


(4*04) 


In  D, 


(2) 


N'(O.k)  =  [ua(0)  0(0)] 1/3  H.j/3[u2(0)] 


f?Q(  0)1  i/? 

2 -  e 

77 


xp 


kv  a( a3) 


2m' 

3 


cos 


3rT| 

kx  (0)  -  — 1 

4 


(4-05) 


In  G, 


N'(0,k) 


f20(  0)  ] 

1/2  f 

l  -  J 

"pn 

TT' 

u(0)  +  — 

12 


} 


(4*06) 


From  Eq.  (4*03)  we  see  that  N'(0,k)  has  no  ieros  for  /3(k)  <  c  ,  so  that  the  series 
(2*10)  begins  with  values  of  kn  such  that  /3< k ^ )  >  c(.  From  Eq.  (4*04)  it  appears 
that  there  are  no  real  roots  in  the  range  Cj  <  |/3(kB)|  <  c0>  It  should  be  noted, 
however,  that  the  asymptotic  expansions  of  the  Bessel  functions  are  valid  only  when  the 
quantities  kx  and  ky  are  large.  Consequently  expf-ky^O)]  is  very  small  compared  to 


■  •wruyr'** 
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exp[kyj(0)]  and  the  values  of  kB  in  this  range  will  be  given  very  nearly  by 

coa  k  x,( a  )  =  0  or 

°  2  1 


W  x  (a  )  -  (n - )  77 

«  a  i  2 


(  4-07) 


whefe  n  is  a  positive  integer.  From  Eq.  (4*05)  the  roots  in  the  range  cQ  <  /3( k )  < 
are  given  by 


*‘n*a'0> 


(n - >77 

4 


(4*08) 


and  from  Eq.  (4*06)  there  are  no  roots  for  /5(  k )  >  c2 . 

Equations  (4*07)  and  (4'08)  for  ,  being  derived  from  the  asymptotic  expansions, 
are  not  valid  for  small  values  of  n,  although  they  are  correct  to  three  significant 
figures  or  better  for  n  -  4.  Explicit  expressions  for  the  f  /wer-order  roots  in  domain  F. 
cannot  be  derived  readily,  but  in  domain  D  they  may  be  computed  from  the  equation 

H-2/3tua(0)1  =  7^  -nl/6  -  J.  i/ 3  [kx  j(  0)  ]  .  (  4’  09) 

The  first  five  roots  of  this  expression  are  kBx2(0)  r  0.6855,  3.9028,  7.0549,  10.2007, 

and  13.3445.  For  frequencies  hif^ier  than  about  0.2  cycle  the  sum  (  Eq.  2*10))  for  cfc  in¬ 
cludes  a  large  number  of  high-order  modes,  and  the  contribution  from  the  few  low-order 
modes  is  comparatively  small.  For  very  much  lower  frequencies  the  only  normal  mode  solu¬ 
tions  that  exist  are  those  of  low  order.  At  such  frequencies,  however,  the  representation 
of  the  function  N(z,k)  in  terms  of  Bessel  functions  becomes  a  poor  approximation.  The 
present  theory  therefore  may  be  considered  valid  only  for  frequencies  sufficiently  high 
for  the  major  contribution  to  <t>  to  come  from  modes  of  order  n  >  4.  For  these  modes 
the  asymptotic  forms  of  the  Bessel  functions  provide  a  sufficiently  accurate  representa¬ 
tion,  and  these  forms  accordingly  will  be  used  in  all  subsequent  development. 


5.  ASYMPTOTIC  APPROXIMATION  FOR  <t 


We  now  have  explicit  expressions  for  all  the  quantities  appearing  in  Eq.  (2*10)  except 
the  factor  [3N'  ( 0 ,  k )  /Bk ]  k-k 

n 

From  the  definition  (3*15), 


_d_ 

dk 


kxa(ai)] 


c  j  <  /S(k)  <  c 0 


(  5*01) 
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For  /P(k)  >  cfll  at  is  replaced  by  0  in  this  expression. 

From  Eqs.  (4*04),  (4*05),  and  (5*01)  with  the  condition  N'(0,k  )  =  0,  we  find 


BN'(O.k) 


Bk 


2k  (-1)“*  1  exp 


k„y2(a3)  +  knVl(n)  " 


rn'i 

12 


(20  f  n ) 


77 


k«k 


f 


1/2  |  d7. 

Q 


_  (5*02) 


for  Cj  <  /3(ko)  <  ce,  and 
BN'(0,k) 


Bk 


2kn(  -1 ) "  exp 


k=k 


■  2771^0.(0) 

kny2(a3)  +  — - 


77 


1/2  I**  d  7 

J.  ^ 


(5-03) 


for  cfl  <  /3(ko)  <  c3.  From  Eq.  (3*20)  combined  with  (5*02)  and  from  (3*10)  with  (5*03)  we 


h  ave 


N(0,k  ) 


rBK,,(o,k)-| 

Bk 


k  =  k 


fik  (  7ri3  |  ‘  dz 

— ^  exP|2kny2(a3)  +  — 


r 


(5*04) 


0. 


for  Cj  <  /3(kn)  <  cQ  and  the  same  expression  with  0  replacing  Bj  for  c0<-  >S( kn )  <  c 

We  shall  assume  the  source  to  be  at  an  altitude  less  than  that  of  the  Telocity  mini¬ 
mum,  that  is.  h  <  i  ,  so  that  N(h,kn)  is  given  by  Eqs.  (3*19)  and  (3*20). 


N(h,k.) 


' 

f  2  I  1/2  ntl 

' 

STril 

' 

2 

(  1 )  exp 

k>n(h)J 

kny2(a3) 

-  k  y,(0)  +  -  cosh 

“  1  12  J 

kn {y j ( h )  -  yj(o)> 

c  <  IHV  )  <  c(h) 

i  n 


tjQ„  ( h ) 


1/2 


exp 


?77i') 

knV2(a3)  +  ~\COa 


f  773 

wh>  -  - 

4j 


c(h)  <  A*(k  )  <  c 

n  2 


(5*05) 


Substituting  Eqs.  (5*04)  and  (5*05)  in  Eq.  (2*10)  and  using  the  asymptotic  form  of 
(  2 ) 

H0  (knr),  we  may  write 

0  =0+0  -  +  <£, 

r  •  •  1  2  j 
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where 

<£,  =  ttA  exp [i(cut  -  n/6]  (“l)"41 

e ><e (  ®  ) 


f..  1 _ 1 1/3  eyp^~iknr  -  k.y2(«3)  ~  k„y,(o)] 


dr/0„ 


coahjk^ty ,(h)  -  y^Oj^N  (z,kn)  ,  (5*06) 


<p  ~  nK  exp[i(oj t ♦  tt/12}] 


1  ^  1/3  exp[~iknr  -  k„y2(a3)] 


e(li)<4(kii)<e( 


(ri..0„(K) 


r; 


d*/Q„ 


cos  [k  x  (h)  —  tt/  4)N(  at , k  ) ] 
0  2  ■ 


(5*07) 


<t>,  -  TrK  exp[i(cJt  +  77/12)] 


V*k-><<:3 


1  -li/J  exp[-ikor  -  kny2(a3)] 


rk.Ojh) 


s: 


dz/o„ 


cos[koxa(h)  -  7T/4]N(x,kn) 


(5*08) 


Now  each  terra  of  <£j  contains  the  factors  exp[~kBy3(a3)]  and  exp  [~kny ,  ( 0)  ]  ,  both 
of  which  are  very  small.  Except  when  [z,/3(kB)]  fails  in  domain  A,  the  factor 
e*p[“k0y2(0]  is  cancelled  by  its  inverse  in  N(z,kB),  but  the  factor  exp[~kBy  j(0) ]  re¬ 
mains  in  all  cases.  Thus  <fil  is  negligibly  small  for  all  values  of  z.  is  likewise 

negligible  because  of  the  factor  exp [~kBy2 ( «3 ) ]  when  [z,  /3(k’B)]  is  in  domain  A;  and 
all  terms  in  4>2  for  which  [z,/^^)]  falls  in  domain  C  are  also  negligible  because  of 
the  factor  exp[-k  y  ( x )]  which  appears  in  N(z,k  )  in  this  domain.  When  [z,/3(k  )]  falls 
in  domain  E,  Eq.  (3*20)  gives 


N(z,k  )  — 


,f  2 


'K  U)J 


1/2 

(-1)  exp 


kBy2(®3)  ~  k„y,(0)  +  cosh  kD{y1(z)  -  y,(0)}j 


and  the  corresponding  terms  of  02  are  negligible  because  of  the  factor  exp[-kny|  (0)]  . 
The  surviving  terms  of  <£Jf  therefore,  are  those  for  which  [z,/3{kn)]  lies  in  domain  D. 
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For  these  we  hare  from  Eq.  (3*15) 


<^2  ~  2A  exp[i(a)t  +  377/4)] 


2n 


*(‘»><£(kD><e0[  " 


-V'1  cos[k  x  (h)  -W/4)X 

rk.O.(h)0.(x)j  J  'dl/<5„ 


cos[knx2(r)  ~  77/4]S[x,a,(kn),aa(kn)] 


(5-09) 


where  S[z,a  (k  ),a  (k  )]  is  a  factor  that  has  the  value  1  lor  a,(k  )  <  z  <  a„(k  ) 

»  o  i  n  \  n  2  n 

and  0  for  z  outside  this  range.  A  similar  argument  leads  to 


~  2A  exp  [i  (oX  +  377/4)] 


C0<^kJ<C2 


2tt 


rkn0n(h)0.(z) 


i/J  exp(iknr)  co3[k  x  (h)  -  -77/4]  x 

■ '  n  2 


r 


dz/0n 


cos  [k  x  (z)  -  77/4]  s[z,0,  a2(kn)] 


(5-10) 


In  Eqs.  (5*09)  and  (5*10)  the  values  of  kn  are  to  be  determined  by  Eqs.  (4*07)  and 
(4*08),  respectively. 


6  .  GEOMETRICAL  INTERPRETATION  OF  THE  NORMAL  MODES 

It  will  be  shown  that  the  summations  indicated  in  Eqs.  (5*09)  and  (5*10)  raj  be 
carried  out  by  replacing  the  sums  by  equivalent  integrals,  which  then  may  be  evaluated 
approximately  by  the  method  of  stationary  phase.  Before  proceeding  with  this  development, 
however,  it  will  be  instructive  to  consider  the  significance  in  terms  of  ray  geometry  of 
the  functions  Q(z,k)  and  x2(z,k)  that  express  the  dependence  of  <$>  on  z.  Let  6(  i.) 
be  the  angle  between  a  given  geometrical  ray  and  the  horizontal.  Then  Snell’s  law  is 
c(z)  sec  d(z)  ~  constant.  Since  the  function  x2(z)  is  defined  only  for  values  of 
Ak)  >  c(  *),  Ak)  is  a  possible  ray  constant  and  defines  a  ray  for  which  sec  (9(z,k)  - 
/S(k)/c(z).  Equation  (3*01)  for  this  case  may  be  written  in  the  form 

Q(z,k)  =  k  v/sec  J$(  z,  k)  “  1  _  k  tan  #{z,k)  •  (6*01) 

The  zeros  of  Q(z,k)  „>re  thus  the  altitudes  at  which  the  corresponding  ray  is  horizontal. 
From  the  second  of  equations  (3*15)  we  have 
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*,(  *) 


r 


tan  #(z,k)  dz 


and  writing 


this  takes  the  fori 


tan  6 


flk) 

c(  z) 


CSC  6  -  cot  8 


~,U)  =  /3(k) 


J'1 


cot  6?(z,k)  dz 


The  first  integral  in  this  expression  is  the  travel  time  a  1  cn g  the  given  ray  from  the  al¬ 
titude  z  to  the  apex  o.  the  ray  at  altitude  ij(  and  the  second  is  the  horizontal 
distance  traversed  (Fig.  4).  Calling  these  qumtities  tj(z,k)  and  rJ(z,k),  respective¬ 
ly,  we  wri te 


Fu  rthe  rmore , 


x3(z,k)  =  /^klt^z.k)  -  r}(z,k) 


cot  6  dz 


ra(  *,>) 


(6*02) 


(6'03) 


The  general  character  of  the  function  *2(z,k)  for  the  type  of  velocity  altitude  function 
assumed  is  shown  in  Fig.  5. 


By  writing  the  factor  cos[kx2(z)  ~  77/4]  in  exponential  form,  it  may  be  seen  that  each 
term  in  Eqs.  (5*09)  and  (5*10)  represents  a  superposition  of  a  pair  of  progressive  waves 
whose  phases  are  given  (apart  froai  a  constant  term)  by 


RAY  HORIZONTAL 


•  mo  i  i  « 


Fi*.  4.  Ray  path  in  (r,x)  plane. 
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Fir.  5.  r'«e«»ral  character  of  the  function. 
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Fig.  6.  Configuration  of  normal  mode  wave  front*. 
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and 


a ) 


t  - 


«k.) 


77 


/Hkm) 


Thus  the  parameter  /3(ka)  ia  to  be  interpreted  as  ths  phase  velocity  of  the  nth  normal 
mode,  and  k^  =  2tt/Xb  where  \  is  the  corresponding  wavelength  (measured  horizontal¬ 
ly).  Tlie  configuration  of  the  wave  fronts,  or  loci  of  constant  phase,  is  illustrated  in 
Fig.  6.  This  configuration  is  such  that  the  velocity  of  propagation  of  any  portion  of  a 
ware  front  measured  along  the  normal  is  c(z),  while  the  velocity  measured  horizontally  is 
the  same  for  all  values  of  z  and  is  equal  to  /3(kQ). 

The  conditions  (4'07)  and  (4*('8)  that  determine  the  k’s  may  be  written  in  the  form 

2xJ(a1)  =  \B(n  - -\)  c(h)  <  /3(ko)  <  cQ  ( 6*04) 

2x3<0)  =  A.(n--T)  c#  <  /S(kB)  <  c,  •  (6’05) 

These  are  the  conditions  which  must  be  satisfied  if  each  of  the  two  sets  of  wave  fronts  is 
to  be  in  phase  with  the  reflected  or  refracted  image  of  the  other,  with  allowance  for  a 
quarter-wavelength  phase  shift  when  the  wave  front  is  reversed  by  refraction  and  no  phase 
shift  when  the  reversal  is  by  reflection  at  the  ground  surface.  This  is  exactly  analogous 
to  the  condition  for  constructive  interference  between  multiple  internal  reflections  that 
determine  the  normal  modes  in  the  case  of  a  medium  made  up  of  discrete  layers.5 

We  also  observe  that  the  (asymptotic)  vanishing  of  <t>  ,  which  corresponds  to  the 
non-excitation  of  the  normal  modes  for  which  /3(kn)  <  c(h),  is  the  expression  in  the 
normal  mode  representation  of  the  fact  that  the  corresponding  rays  do  not  pass  through 
the  level  of  the  source.  The  interpretation  of  /3(k),  for  real  values  of  k,  as  a 
quantity  defining  a  family  of  rays  in  the  (r,z)-plane  indicates  the  physical  meaning  «f 
the  integrals  around  the  cuts  in  the  k-plane.  The  integral  around  the  cut  starting  at 
cj/c  arose  by  transformation  of  the  original  path  of  integration  along  the  real  axis 
for  0  <  /3(k )  ^  c  .  The  corresponding  rays  are  confined  to  the  portion  of  the  (r,z)-plane 
for  which  z  is  greater  than  the  altitude  a3  (Fig.  7).  Such  rays,  therefore,  do  not 
contribute  appreciably  to  the  value  of  4>  when  either  h  or  z  is  less  than  aJ(  which 
is  the  case  in  which  we  are  primarily  interested.  Similarly,  the  integral  around  the  cut 
starting  at  ^/c2  represents  the  contribution  to  <t>  of  rays  for  which  /3(k)  >  c  , 
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Fig.  7.  Rays  for  which  0  <  /Vk)  <  c. 

Fir.  R.  Ri  ys  for  which  . 

that  is,  of  the  type  shown  in  Fig.  8.  These  rays,  after  once  passing  above  the  level  z^, 
thereafter  do  not  return  to  lower  elevations  and  thus  do  not  contribute  appreciably  to  <t> 
when  the  point  (r,z)  lies  to  the  right  of  the  upward  refracted  branch  of  the  critical  ray 
that  becomes  horizontal  at  the  altitude  z^.  For  such  points,  therefore,  repre¬ 

sents  essentially  the  whole  solution. 


7.  PHASE  AND  GROUP  VELOCITIES  OF  THE  NORMAL  MODES 


Since,  for  a  given  value  of  n,  /3(kR)  *s  a  function  of  frequency  defined  by  F.qs. 
(4*07)  or  (4*08),  the  individual  terms  of  Eqs.  ( 5  *  09 )  and  (5*10)  represent  disperseo 
waves.  Therefore,  it  is  of  interest  to  compute  the  corresponding  group  velocities.  The 
group  velocity  Un  corresponding  to  the  phase  velocity  -  /3(kn)  is  defined  by 


U  =  da;  /dk 

o  n 

From  Eq.  (4*07)  we  have  for  n  =  constant 


(7-01) 


X2(ai)dk„  +  kn  dx3(a1)  ~  0 

or 


X2(8|) 


Bx2(  )' 

■ 

3k 

n  J 

w 

dk 


f3x2(a1)T 


d  u> 


O 


From  the  fact  that 


(7*02) 


x 


-  k2  dz 

II 


mmm  m*  4 


2ft 


and  the  fact  that  the  integrand  vanishes  at  » (  and  aJ(  we  have 


rBx2(ai)i 


3o> 


and 


rf  - 


dz 


(z)  /  1  - 


c(  z) 


l/^n) 


r 


dz 

c(  z)sin  (9 


tJUi)  ( 7  ’03) 


3X2(aI) 


3k 


whence 


r 


—  - 1  tan  6  dz  —  - 1  cot  0  dz 

k  k 


a: 


X2(ai}  +  r2(8, 


(7-04) 


_  ^nS(ai) 


l) 


r7(Bl) 

s<v 


(7-05) 


Equation  (4*08)  leads  to  the  same  relations  with  aj  replaced  by  0.  Thus  the  group 
velocity  of  the  nth  mode  is  equal  to  the  average  horizontal  velocity  along  the  corres¬ 
ponding  ray  between  its  maximum  and  minimum  altitudes. 

A  numerical  example  is  plotted  in  Figs.  9(a)  and  9(b)  which  show  /3n  and  Un  as 
functions  of  frequency  for  the  following  velocity-altitude  function: 


Alt.  (km) 


c ( m/ sec ) 


0 

16.  5 

32 

50 


344 

286 

3)2 

380 


Linear  velocity  gradients  are  assumed  between  the  tabulated  values,  arid  the  velocity  is 
assumed  to  decrease  monotoni ca 1 ly  above  50  km.  The  fact  that  this  velocity  function 
does  not  fulfill  the  condition  of  continuity  of  dc/dz  is  beside  the  point  since  the 
ray  geometry  would  not  be  greatly  altered  i  f  we  supposed  the  discontinuities  in  gradient 
to  be  rounded  off  to  satisfy  the  continuity  condition.  The  phase  velocity  curves  provide 
the  basis  for  the  statement  previously  made  that,  for  frequencies  higher  than  about  0.2 
cycle,  c i)  includes  a  large  nin&ber  of  high-order  modes.  For  example,  the  number  of  terms 
in  <ti  for  any  given  frequency  is  ecm»l  to  the  number  of  values  of  Bn  at  that  frequency 
falling  between  the  limits  cQ  anti  Thus  at  0.2  cycle  there  are  12  terms,  from  the 
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22nd  to  the  33rd  mode,  inclusive;  at  1.0  cycle  there  are  58  terns,  from  the  107th  to  the 
164th  mode;  and  at  5.0  cycles  there  are  288  terms,  from  the  532nd  to  the  819th  mode. 


8.  INTEGRAL  REPRESENTATION  OF  <t> 


Since  the  number  of  terms  is  so  large,  Eqs.  (5*09)  and  (5*10)  are  not  very  useful  as 
they  stand.  However,  the  fact  that  successive  values  of  kn  differ  only  slightly  sug¬ 
gests  that  the  sums  may  be  converted  into  equivalent  integrals.  Using  the  symbol  A  to 
denote  the  change  in  any  quantity  from  the  nth  to  the  (n  +  l)th  mode,  we  have  from 
either  F.q.  (4*07)  or  Eq.  (4*08) 


[1/7t]  a  [kn*2(a1>]  -  An  -  1 


Bu  t 


A[k„xa<ai>]  '  Akn{x,<a,>  +  k  n  — 3k->  } 


-  r2( a j )  Akn 


from  F.q.  (7*04).  Hence 


~ra(a.) 


A  k 
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i8*01) 


Multiplying  each  term  of  Eqs.  (5*09)  and  ( 5  *  10 )  by  (8*01)  and  using  Eqs.  (6*01)  and  (6*03), 
we  have 


<t>2  -  2A  exp  [ 


i(cut  -  77/4] 


C<h><^n<e0 


77rk  tan#  (  z)  tan  ^(h) 

n  •  ' 


l  1/2  exp(~iknr)  cos[kBx2(h)  -  77/4]  x 


cos[knx2(z)  ~  77/4]  S(z,  a  j ,  a2)Akn  (8*02) 


am 


<t>  ~  2A  exp[i(c*jt  77/4] 


e  <0  <e, 

0  n  2 


[Trrk|itan  8Q(z)tan  #n(h)J 


1/  3 


exp(-iknr)  cos[knx2(h)  -77/4]* 


cos[k0x2(z)  -  tt/4]  S(  z,0,  a2)Akn 


(8*03) 
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Equation  (8*01)  shows  that  when  rj(aj)  is  a  lflrRe  number  of  wavelengths,  Akn  may  be 
treated  as  a  small  quantity.  However,  before  passing  to  the  limit  Ak^  — •  dkn  and 
writing  these  sums  in  the  form  of  integrals  over  a  continuous  range  of  values  of  ko,  we 
must  consider  that  in  deriving  Eqs.  (5*09)  and  (5*10)  terms  depending  upon  x2(at)  and 
x  ^ ( 0 )  were  eliminated  by  use  of  the  relations  ( 4  *  07 )  and  (4*08).  Refore  allowing  kn 
to  assume  continuous  values  we  reintroduce  the  dependence  upon  x2(a])  and  x^O)  by 
noting  that  for  any  integer,  s, 

exp{-2is  [knx2(a  j )  +  tt/2]}  =  exp(~2i77ns)  =  1  c(h)  <  /3m  <  cQ 

and 

exp{-2is[k(|xJ(0)  +  377/4]}  =  exp(~2i77ns)  =  1  co  <  <  ca 


These  expressions,  therefore,  may  be  introduced  as  factors  after  the  summation  signs  in 
Eqs.  (8*02)  and  (8*03).  The  reason  for  the  introduction  of  the  arbitrary  integer  s  will 
become  apparent  presently.  Passing  now  to  the  limit  A  kn  — *  dk  and  writing  the  cosine 
terms  in  exponential  form,  Eqs.  (8*02)  and  (8*03)  become 

<t>3  =  —  exp[i(oX  -  tt/4  -  stt)]  (e'"*72^  +  e™72^  +  I3  IJ  (8*04) 

an  d 

<PS  -  —  exp[i(oJt  ~  77/4  “  3377/2)]  (e_TTi/  2  I5  +  e77^2  +  I?  +  Ig)  .  (8*05) 

whe  re 


♦ 

0 

f  ( k )  exp{-  ik  [  r  -  x^h)  -  x2(z)  +  2sxa(a.)]}  S(  * . «,  *  )  A 


(8*06) 


f(k) 


exp{~  ik  [r  +  x  2(  h )  +  Xj(z)  +  2sxa(a,)]}  5(z,a1,aj)  dk 


(8*07) 
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I 


s 


I 


4 


exp{-  ik[r  -  xf(h)  +  x}(i)  +  2axJ(»1)]}  ?(*,*,. *a)  A 


exp{-  ik [r  +  x2(h)  -  x2(z)  +  2sxj(«i)]}  Sfz.a^aj)  dk 


'8  *08) 


(8*09) 


f(k) 


2 


'I  1/  2 


[rrrk  tar  #(z,k)tan  £(h,k) 


(8*10) 


The  lower  limit  of  in  tegration ,  kh,  is  the  first  value  of  hn  less  than  to/ c(h) ,  and 
k*  is  the  value  of  ka  next  larger  than  ^/cQ.  The  integrals  I?  to  I#  are  given  by 
the  same  expressions  with  aj  =  0  and  the  lower  and  upper  limits  of  integration  re¬ 
placed  by  k"  and  kjf  respectively,  where  k“  is  the  value  of  kn  next  smaller  than 

oj/c  and  k  „  is  the  value  next  larger  than  a'/c  . 

'02  n  2 

9.  APPROXIMATE  EVALUATION  OF  THE  INTEGRALS 


In  the  approximate  evaluation  of  integrals  of  the  type  (8*06)  to  (8*09),  we  must 
consider  the  criteria  for  the  existence  of  stationary  points  and  points  of  inflection  of 
the  phase  of  the  exponential  factors  in  the  integrands.  We  let  X(z,h,k)  stand  for 
any  one  of  the  four  quantities  2axJ(a1)  ±  x2(z)  ±  x2(h),  so  that  the  phase  may  be  writ¬ 
ten  in  the  form 

( r,  x,  h,  k)  =  k[r  +  X(  z,h,  k)]  (9*01) 

We  also  introduce  the  notation  R(z,h,k)  =  2sr^(*l)  ±  r2(z)  ±  r2(h)  and  T(r,h,k) 
=  2st2(aj)  ±  t2(z)  ±  t  a  ( h ) .  Then,  using  Eqs.  (5*01)  and  (6*03)  (which  remain  true  if  i 
or  h  is  substituted  for  the  lower  limit  of  integration,  a  )  we  have 
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<3  y; 

3  k 

r  -  R(z,h,k) 

(0-02) 

B20 

BH 

(0-03) 

-  z: 

—  - 

Bk2 

Bk 

B30 

B2R 

(9-04) 

B  k3 

Bk2 

The  condition  for  a  point  of  stationary  phase  is  then  r  ~  R(z,h,k),  which  is  to  be  re¬ 
garded  as  an  equation  to  determine  a  particular  value  (or  set  of  values)  of  k  for 
given  values  of  r,  z,  and  h.  The  geometrical  significance  of  this  condition  is  illus¬ 
trated  in  Fig.  10,  which  shows  a  number  of  refracted  cycles  of  the  pair  of  rays  specified 

by  a  particular  value  of  k,  say  K.  The  radial  distances  to  the  points  Aj  ,  A2 . 

B  ,  Bj  ,  etc.  at  altitude  z  are  given  by  the  equations, 


"  A 

2sr 2( a  j ) 

+  r  j(  h ) 

1 

N 

s 

r  n, 

i. 

2 . 

=  2srj(ai) 

+  r j(h ) 

+  r2(z) 

s 

=  o, 

i. 

2 . 

rc 

=  2srj(aj) 

~  r2(h) 

1 

K) 

N 

s 

=  l. 

2, 

3 . 

r  D 

=  2sra( a  j ) 

-  r2(h) 

+  r2(z) 

s 

=  l. 

2, 

3 . 

all  of  which  are  comprised  in  the  general  condition  for  the  existence  of  a  point  of  sta¬ 
tionary  phase  for  one  or  more  of  the  integrals  (8’06)  to  (8*09).  Thus  we  see  that  when¬ 
ever  the  point  of  observation  (r,z)  lies  on  a  possible  geometrical  ray  from  the  source, 


Fig.  10.  Ray  geometry  and  points  of  stationary  phase  . 
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at  least  one  of  the  integrals  Ij  to  I#  will  have  a  point  of  stationary  phase  for  some 
value  of  s.  If  (r,z)  falls  at  a  point  where  two  rays  heving  different  values  of  K  in¬ 
tersect,  there  will  be  two  points  of  stationary  phase;  and  if  (r,z)  lies  in  a  geometrical 
shadow  zone,  the  equation  r  -  R(z,h,k)  will  have  no  solution  for  real  values  of  k  and 
none  of  the  integrals  will  have  points  of  stationary  phase.  The  integer  s  is  now  seen 
to  represent  the  number  of  times  the  given  ray  passes  through  its  minimum  altitude  (or  is 
reflected  at  the  ground  surface)  in  reaching  the  point  (r,z). 

Approximate  evaluations*  of  the  general  integral  appropriate  under  various  circum¬ 
stances  are  as  follows: 


9.  1  CASE  I 

(r,z)  does  not  lie  in  a  shadow  zone  and  is  not  a  focal  point,  so  that  r  -  R(K),  but 
3R/3K  f  0. 


f-~] 

3R/3K) 


1/3 


f(K)  exp{-iK[r  +  X(K)]}  I  exp  [~i77va/2]  dv 


J. 


(9-06) 


whe  re 


+ 


£ 


K[r  +  X(K)  ]  -  k.[r  +  X(k.)] 


1/2 


(9-07) 


'2  r  ,11 »/» 

far  K[r  +  X(K)]  -  ka[r  +  X(k.)] 


(9*08) 


In  the  integrals  Ij  to  I4,  kfe  -  k*,  k,  -  kh,  and  in  the  integrals  Is  to 

T  lc  —  W  L  ~  k 

ls>  *b  *2’  •  *0# 

If  the  approximate  evaluation  of  I  represented  by  Eq.  (9*0&)  is  carried  further, 
the  next  term  is 


Jf(kb)  -  f(K) 
Lr  -  R(kb) 


exp{~ikb[r  +  X( kfc ) ] } 


f(k.)  -  f(K) 

r  -  R(k.) 


exp{-ik.[r  +  XU,)]}} 


8S»*.  for  •■■■pic,  t  ha  treatment  of  tha  oat  hod  of  atationary  pbaaa  |lran  by  C-  Eckart,  "Approiiaata 
Solution*  of  Ona  Dlatm  tonal  »a.a  Iquatlona."  Par.  Hod.  pbya..  rol  20,  p  399:  1948- 
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This  term  should  be  small  if  Eq.  (9*06)  is  to  be  used  aa  it  stands.  For  a  given  value  of 
/3(k),  f(k)  and  BR/Bk  are  proportional  to  so  that,  apart  from  the  fluctuations 

due  to  the  presence  of  the  Fresnel  integral,  the  absolute  value  of  the  first  term  is 
independent  of  the  frequency  and  the  second  is  proportional  to  In  addition  to 

this  explicit  dependence  on  the  frequency,  however,  there  is  an  implicit  dependence  due 
to  the  fact  that,  with  increasing  frequency,  kh  approaches  cu/c(h)  so  that  tan  #(h,kh) 
approaches  zero  and  f(k)h  becomes  very  large.  Therefore,  we  have  no  assurance  that  the 
second  term  will  be  small  for  any  frequency.  This  convergence  difficulty  arises  from  the 
fact  that  the  asymptotic  approximations  used  for  the  function  N(z,k)  are  not  valid  for 
values  of  z  and  k  such  that  tan  6  is  very  small.  If  it  were  possible  to  obtain  exact 
solutions  of  the  differential  equation,  the  convergence  difficulty  presumably  would  not 
occur.  Therefore,  we  take  Eq.  (9*06)  as  it  itands,  without  any  correction  term,  as  repre¬ 
senting  the  best  approximation  obtainable  by  the  present  method.  As  justification  for 
this  pro»  edure  we  shall  show  that,  except  for  the  diffraction  effects  represented  by  the 
Fresnel  integral,  Eq.  ( 9*06)  leads  to  exactly  the  same  expression  for  the  energy  flux 
density  that  one  derives  on  the  basis  of  ray  geometry. 

When  vfc  and  vf  are  large  in  absolute  value,  the  Fresnel  integral  approaches  the 
value  (2i)1/'  so  that 


I 


2 


r.  br 

i  r  K  — —  tan  8  ( z ,  K )  tan 


J/  2 


exp{-iK[r  +  X(K)]} 


(9*09) 


Calling  <t>  (K)  the  contribution  to  that  arises  from  any  particular  ray  through  the 

point  (r,z),  we  have 


Amp  <t>  (K) 


B  R 

RK - tan  0  (? ,  K)  tan0(h,  K) 

B  K 


1/2 


(9'  10) 


In  computing  the  particle  velocity  from  we  may  treat  Amp  <t>  as  constant  and  differen¬ 
tiate  only  the  exponential  phase  factor,  since  Amp  <f>  u,  by  conqparison,  a  slowly  varying 
quantity.  By  use  of  Eq. (1*07),  the  absolute  value  of  the  energy  flux  density  is  then 


|J(K) 


P(  h) 


c(  z) 
p(h)co3  A3 


[go  )]  2 

BR 


c(z) 


RK 


BK 


tan  8  (  z,  K)  tan  8  (h  ,  K) 


-  l 


(9-11) 


Fro*  the  viewpoint  of  ray  geometry,  we  consider®  the  energy  flux  emitted  in  the  aolid 
angle  between  the  rays  that  leave  the  source  at  angles  8  ( h,  K)  and  0  (h.K)  +  68  (Fig.  11). 
If  E  is  the  total  rate  of  energy  emission  at  the  source,  we  have 

F  F 

d  E  -  —  cos<9(h,K)d$  -  —  d[sini9(h,K)] 

2  2 

The  area  of  wave  front  ircluded  between  these  rays  at  the  point  (z,R(K)]  is 


dS 


~  2  77  R  sin  8  (z,K)  dR 


Hence 


dE 

dS 


r  dR 

F.  I  JJ—  0(-h— j 


t9* 12) 


Nov 


dR 


d  [sin  6  (h,  K)  ] 


BR 

BK 


B  [sin  8  ( h  ,  K ) ) 

BK 


and  using 


sin  6  (h,  K) 


1  -  KJc3(h)/  to  2 


1/2 


with 


c  os  8  (z,  K)/cos  8(  h,K)  -  c(z)/c(h) 


Fa.  (9*12)  may  be  written  in  the  form 
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j  RK  — —  tan#(z,K)  tan  0  ( h.K) 
4 7 T  c(  z)  [  B  K 


(9-13) 


From  the  limiting  form  approached  by  (fi  in 
the  neighborhood  of  the  point  source,  F 
and  A  are  connected  by  the  relation 


477/?  (h)  co  2  A2 
c(h) 


(9-14) 


which  establishes  the  identity  of  Eqs. 
(9*13)  and  (9*11).  Equation  (9*06)  thus  has 
been  shown  to  be  completely  equivalent  to 
geometric  ray  theory  in  the  high-frequency 


Thia  (rtatatnt  folio**  a  aiailar  derivation  by  B.  Gatanbarf,  in  "Propagation  of  Sound  Vara*  In  tha 
Atwoaphera,"  J  A  ■"  o  □  *  t  So  c  .Am  ,  vol  14.  p  1 S  1 .  1942-  It  1*  rapaatad  hare  tinea  the  identity  be¬ 

tween  lutenberi'i  equation  and  tq.  (9*11)  i*  »  >t  oba  lout  without  tone  change*  in  notation. 
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The  above  argument  also  shows  th*t  Eqs.  (1*05)  and  (1*06)  conform  to  the  require¬ 
ments  of  energy  conservation. 

9.2  CASE  II 

(r,z)  lies  in  a  shadow  zone  where  the  boundary  of  the  shadow  zone  is  formed  by  a  single 
ray  (i.e.,  the  boundary  is  not  a  caustic). 

Let  K  be  the  parameter  specifying  the  bounding  ray.  Then 


I  -  f(K)  I  exp{-i[a  +  b(k  -  K)  +  c(k-K)J]}dk 


rhere  a  -  K[r  +  X(K)],  b  -  r  ~R(K),  c  -  -%  BR/BK.  Substitution  of 

i  =  (2/ttc)  1/3  [c(k  -  K)  +  b/2] 

ts  the  variable  of  integration  leads  to 


I 


77 


br/bkJ 


i/ a 


f(K) 


exp 


K[r  +  X(K) ]  + 


[r  -  R(K)] 
2  BR/BK 


exp(  ~i  77  va/2) 


dv 


(9-15) 


where 


r  ~  R(K)  ~  (k,  -  K)  BR/  BK 

/-7t  BR/BK  *  (9-16) 


and 


r  -  R(K)  -  (kb  -  K)  BR/BK 
BR/BK 


(9-17) 


9.3  CASE  III 

(r,z)  lies  near  a  point  such  that  BR/ dk  :  0  is  satisfied  for  some  value  of  k  =  K. 

The  appropriate  expression  for  I  in  this  case  is 


;  6 
dJR/dKJ 


1/3 


f(K) 


exp{~iK[r  +  X(K)] } 


b 


I 


exp[i(wv  +  Vs)]  dv 


(9-18) 
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wh  ere 


[R(K)  -  r] 


6 

92R/BK2. 


1/3 


(9-19) 


V 


b 


-32R/BK2-i 

6 


i/ j 


(k.  -  K) 


(32  R  /  3 K 2'i  i/3 

HH  (k‘ ' K) 


(9-20) 


(9-21) 


When  vb  and  ~  y m  are  large  in  absolute  magnitude,  the  integral  in  Eq.  (9*18)  approaches 
the  Airy  integral  (except  for  a  constant  factor), 


When  (r,z)  lies  in  a  shadow  zone,  w  >  0  and  the  Airy  integral  decreases  monotonically 
with  increasing  |wj,  When  (r,z)  is  in  a  geometrical  ware  zone,  w  <  0  and  the  integral 
oscillates  with  slowly  decreasing  amplitude  with  increasing  |w|.  When  w  =  0  the  inte¬ 
gral  has  the  value  2^/3T(2/3)  ~  1.5466. 


10.  DEPENDENCE  OF  DIFFRACTION  EFFECTS  ON  RANGE  AND  FREQUENCY 


To  illustrate  the  way  in  which  the  geoatetrical  ly  determined  amplitude  is  modified  by 
diffraction  we  consider  the  case  in  which  the  source  and  point  of  reception  are  at  the 
ground  surface,  r  =  h  =  0.  For  this  case  ~  0  and  in 


I 


3 


exp{-ik[r  +  2(s  “  1  )  Xj(0)J  }  dk 


22) 


( 10*011 
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f  2 

|  f(k)  exp{— 1 k  [  r  +  2(s  +  1  )  x^O)]}  dk  ,  (  :Q.02) 


I7=  !g  =  f(k)  exP{-ik[r  +  2s  x  j(  0)]  }  dk  .  (10*03) 

These  integrals  have  points  of  stationary  phase  for  r  r  2(s  -  1)  r2(0),  2(s  ♦  1)  r^lO), 


and  2sr0(0),  respectively.  For  s  =  0,  only 

has  a  point  < 

}f  stationary  phase;  hut  the 

same  value  of  r  is  a  point  of  stationary  phase 

with  s  -  1.  Therefore,  since 

for  I5  with 

s  =  2,  and  for  I?  and  I# 

e-  (  1  )  e-TT*/  J 

*5.  1 

■-  1  )  ^i/7  » 

1  6  .  • 

-  *(JttI./2)  i 

-  1  C  *7.  • 

Eq.  (8 -05)  becomes 

<t>  -  2A  exp  [l  (cu  t  - 

77/4  -  3s7*/2)]  I 

(10*04) 

where  I  is  given  by 

Eq.  (10*03)  snd  s  = 

1,2,3 . 

Further,  R  r  2  s  r  2  (  0  )  and 

T  =  2stj(0).  Values  of 

2Xj(0),  2ra(0),  amd  2Br2(0)/d/3  for  various  values  of  /J,  cocs- 

puted  for  the  velocity 

-altitude  function  given  above,  are 

TABLE  1 

Values  of  Parameters  Occur r ing  in  I 

given  in  Table  1. 

0 

2xa(0) 

2r  (0) 

(wn) 

2dra(0)/d/3 

(km/sec) 

(km) 

( sec ) 

0.344 

36.562 

240.63 

-  00 

.34449 

37.154 

230.79 

-9088 

.345 

37.344 

227.29 

-5882 

.347 

38.852 

219.51 

-2740 

.350 

41.070 

213.49 

-1494 

.355 

44.668 

208.25 

-720 

.360 

48.208 

205.65 

-358.2 

.365 

51.710 

204. 43 

-144.6 

.370 

55.224 

204.08 

-2.2 

.375 

58.740 

204.  33 

+99.2 

.380 

62.248 

205.03 

+174.4 
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The  infinity  in  dra(0)/d fi  at  /I  r  c  ,  which  is  of  an  order  such  that  tan  $(0,k) 
dr2(0)/d/3  is  finite,  would  lead  to  difficulties  for  r  near  R(kQ)  i  f  we  took  the  lower 
limit  of  integration  in  Eq.  (10*03)  to  be  ^  C,)/C0  instead  of  k^ .  Fxcep*  where  this 

infinity  is  involved,  however,  we  may  take  the  integration  limits  to  be  m/c^  and  w/cj( 
since  these  differ  very  slightl>  from  k^  and  k  .  (For  example,  for  a  frequency  of 
1  cycle,  Akg)  r  0.34449  and  /3( k 2 )  =  0.37961  in  place  of  0.344  and  0.380,  respec¬ 
tively.  ) 

From  Table  1  we  note  that  r^O)  has  a  minimum  for  /3  -  0.370  (nearly).  Case  III, 
with  s  -  1,  2,  etc.,  is  therefore  the  appropriate  approximation  for  values  of  r  near 
204,  408,  etc.  At  these  minima  of  R  we  have  BsR/3/92  -  2.4  x  10*s  se^km*1,  from 

which  we  find  the  l\mits  of  the  integral  in  Eq.  (9*18)  are  vfc  -  “0.  553(  va  )*  /  J  and 

v^  ~  1 . 589(  ^s) 1 /3.  Except  for  quite  high  frequencies,  or  larg  "ilues  of  s,  these  quan¬ 
tities  are  not  sufficiently  large  for  the  Airy  integral  with  infinite  limits  to  be  ap¬ 

plicable.  The  corresponding  value  of  the  parameter  w  is  0.808  (204.08s  ~  r)j^2/3  s‘  1/3, 
and  when  |204.08s  -  r|  >  5.69s,  |w|  >  |v  -  vb|3.  Under  these  conditions  we  can  obtain  a 

crudely  approximate  evaluation  of  the  integral  in  Eq.  (9*18)  by  neglecting  the  v3  term 
in  the  exponent.  The  absolute  value  of  this  integral  is  then 


(2 [ 1  ~  cos  w( va  -  vfa ) ] > 1 7  2 


1  75s1' 3  v  ' 3/3  r 

— - Jl  -  cos  1.73  204.08s 

204.08s  -  r|  \ 


from  which  the  distance  Ar  between  diffraction  maxima  is  3.632A'  km  and  the  envelope 
of  the  maxima  falls  off  as  2.  476s  1  /3/ |  204.  08s  ~  r\u3/i. 

Case  I  is  the  appropriate  approximation  for  the  integral  I  in  the  ranges 
205  +  <  r  <  2ra(0,k^)  with  s  =  1,  410  +  <  r  <  4ra(0,ke)  with  s  =  2,  etc.  Values 
of  v^  v'l/2  and  vfa  u't/3  for  this  case  are  given  for  various  values  of  r  in  Table  2- 
From  these  figures  it  appears  that  for  frequencies  of  the  order  of  1  cycle  the  limits  of 
the  Fresnel  integral  in  Eq.  (9*06)  are  not  large,  so  that  diffraction  minima  and  maxima 
will  occur  over  the  whole  width  of  the  geometrical  wave  zone. 

Case  II  is  applicable  for  values  of  r  greater  than  2sra(0,k‘),  which  for  a  1-cycle 
wave  has  the  values  230.79,  461-58,  etc.  in  the  present  example.  The  limits  of  the  Fresnel 
integral  as  given  by  Eqs.  (9*16)  and  (9*17)  are  both  of  the  same  sign,  so  that  the  integral 
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decreases  aoDOtonical ly  as  r  increases.  Inserting  numerical  values  for  v  1,  Kqs. 

(9’16)  and  (9*17)  become 

r  -  230.79s 
23 . 20s 1/2 

r  +  91 . 6s 
b  23. 20s,/J 

The  upper  limit,  vb,  is  thus  a  fairly  large  quantity  for  all  values  of  s,  and  the  value 
of  the  Fresnel  integral  is  determined  chiefly  by  v^. 


Values  of  Limits  of  Fre 
r(km) 

205-65 

208.25 

213.49 
219.51 

227.29 
240.63 

411.30 

416.50 
426.98 
439.02 
454.58 

481.26 


TABLF.  2 

l  Integral  in  Fq .  (9 ‘06) 

vbA',/J 

0.82 

1.55 

2.57 

3.51 

4.51 
5.91 

1.16 

2.19 

3.63 

4.96 

6.38 

8.36 


for  Var  ious  Values  of  r 
v  /v,/J 

m' 

-2.06 

-1.77 

-1.29 

-0.81 

-0.40 

0.0 

-2.92 

“2.50 

-1.83 

-1.14 

-0.57 

0.0. 


The  resulting  values  of  |d>|  in  and  near  the  first  geometrical  wave  zone  (s  =  1)  are 
plotted  as  a  function  of  r  for  a  frequency  of  1  cycle  (Fig.  12).  The  dashed  curve  shows 
the  •geometrical*  value,  i.e.,  the  limit  approached  for  infinite  frequency.  At  longer 
ranges,  correspond ing  to  larger  values  of  a,  the  values  of  v^  and  vfa  given  by  Eqs. 
(9*07)  and  (9*08)  increase  in  proportion  to  s,/J,  so  that  the  oscillations  of  the  Fresnel 
integral  within  the  geometrical  wave  tone  become  smaller.  The  limits  of  integration  given 
by  Ekja.  (9*20)  and  (9*21)  increase  in  proportion  to  s1/3,  while  the  parameter  w,  for  a 


Fig.  12.  Relative  amplitude  distribution  in  and  near  firat  geometric  wave  zone. 

given  distance  [R  (K)  -  r]  from  the  edge  of  the  geometrical  shadow  zone,  decreases  as  s' 1/3. 
The  complete  Airy  integral  thus  becomes  a  better  approximation  for  the  integral  with 
finite  limits  in  Eq.  ( 9  *  18 ) .  Igioring  the  minor  diffraction  oscillations,  the  gross  vari¬ 
ations  of  \<b\  with  distance  at  very  long  ranges  are  shown  in  Fig.  13,  for  a  frequency  of 
1  cycle.  The  dashed  curve  shows  the  value  of  \<t>\  for  a  source  of  the  same  strength  in  a 
medium  of  constant  propagation  velocity. 


OASmES  Cufcvt  GIVES  AM^liTjOC  VS 
distance  eo«  a  point  solace  o#  tm£ 


Fig*  13.  Relative  amplitude  at  long  ranges. 


